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CROP IDENTIFICATION FROM RADAR I\V\GERY OF THE 
HUNTINGTON COUNTY, INDIANA TEST SITE 


ABSTRACT 


The resulH of a s^uc^y fo discriminafe crop f-ypes using L~band, dual polarizal’ion 
(HH and HV) radar data are reported. X"band data unfortunately were not available 
for analysis due ro problems encountert'd during the flight. The flight was made over 
Huntington County, Indiana on SeptemL-»r 13, 1973 using the ERIM radar. The test 
site consisted of fields of corn, soybeans, voods and pasture. 

The analysis resulted in the following observations: 
a) Like polarization was successful in discriminating corn and 
soybeans, however pasture and woods were consistantly 
confused as soybeans and corn, respectively. The probability 
of correct classification was about 65%. 


b) The cross polarization component (highest for woods and lowest 
for posture) helped in separating the woods from corn,and pasture 
from soybeans and when used with the like polarization component, 
the probability of coTact classification increased to 74%, 



1.0 INTRODUCTION 




* 

On Sepfember 13, 1973 the ERIM synthetic aperture radar was flown over 
Huntington County, Indiana. The test site wr; covered by two passes (Pass 2 and 
Pass 3). Both passes were flown at 6500' msl and had a display width of 19,000' . 

Both X (3.2 cm) and L~band (23.0 cm) imagery for HH and HV polarizations were 
obtained, however the two channels of the X-band data were not supplied by ERIM. 
Two problems , which have been documented in greater detail [1] in ar other report, 
encountered during the flight were responsible for degradation of the X“band image. 
The first problem was the "failure of the interface between the aircraft interial 
navigation system (INS) and the two signal film drives. . . . The second problem 
encountered . . . (was) that an X~band antenna wander problem had occurred ."[ 1] 

The objective of this report is to document the data processing and analysis 
of the L“band radar imagery and report the results of a discriminant analysis performed 
on the data. 


2.0 TEST SITE 


Map 1 shows the areas imaged for both Pass 2 and Pass 3. An aerial photograph 
of the entire test site along with an ASCS overlay indicating field boundaries and 
field numbers were used to obtain coordinates. Figure 1 shows a portion of the aerial 
photcgiaph and overlay . About 70~75% of the fields were corn and soybeans, the 
other two vegetation types (about 2CP/o) were woods and pasture, 5% of the fields 
were unidentified. There was an equal percentage of com and soybean fields . The 
ground truth data supplied were the field numbers, corresponding crop types, percent 
cover and row direction (if applicable). 

The natural vegetation for the area is predominantly beach-maple forest, and 
some oak hickory forest. Both corn and soybeans in September are mature and ready 
for harvest. 


★ 

Environmental Research Institute of Michigan. 
Agricultural Stabilization and Conservation Service. 




40° 41' N 
85° 33’ W 


40° 33.3' N 
85° 29.5' W 


40° 41' N 
85° 29. 8' W 


^ See Figure 1 for Aerial Photograph of Section 


Map 1 . Hunfingtori County test site. 






Aerial photograph and overlay showing section of test site 
1 statute mile . 











3.0 DATA PROCESSING 


Two mognefic fopes wifh radar data of Pass 2 and Poss 3 were mailed to the 
University of Kansas from ERIM on January 17, 1975. The tapes contained only two 
channels of data — l-(HH) and L(HV) — end were registered and interleaved in the 
I-ARSYS III format. The tapes were gener-ted on ERIM's PDP-11/45. The first set of 
rapes sent by ERIM for the Phoenix experiment [2] were generated on the Michigan 
Computing System (MTS). The computational center facilities at the University of 
Kansas (KUCC) encountered unacceptable number of parity errors. A possible explana- 
tion was an incompatibility of the two systems (MTS and KUCC). It was then decided 
that all data to be delivered by ERIM to the University cf Kansas were to be generated 
on ERIM's PDP"ll/45 computer which KUCC had no problem in reading. 

As the HW 635 is a 36-bit/word machine, translitertJion programs were written 
to generate HW 635 compatable tapes [3j . A digital printout of the entire test site was 
produced on which field boundaries were marked. The coordinates of the fields were 
input into the computer, which extracted data from such field and stored the informa- 
tion on tape to be later corrected. 

For an imaging radar, the backscattered return ftom a given type of target 
(such as corn fields) can vary across the image between the near range and far range 
because of the following reasons: (a) antenna gain variations as a function of look 
angle, (b) path loss variation as a function of range and (c) scattering coefficient 
variations as a function of angle of incidence. 

The relative radar response curves (taken from an ERIM report [1]) shown in 
Figure 2 correct for antenna gain and range variation as a function of angle of incidence 
but were not used by ERIM to correct the data due to lack of confidence In them*. To 
inspect the trend of the data the test site was broken up into five strips (parcllel to the 
flight direction of the aircraft) and from each strip all com and soybeans data were 
averaged separately and plotted. The data are shown in Figure 3. The shapes of the 
curves in Figures 2 and 3 are very similar. Corn gave a higher rerum than soybeans 
consistently over tl. -ntire range. 

To correct the data, th • test site was broken up into 250 strips and all soybeans 
data (having me largest number of points) within each strip were averaged and plotted 

*Personal communication, D. Ausherman, ERIM. 
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FIGURE 7. RELATIVE PvADAR RESPONSE, HUNTINGTON, INDIANA, 
PASSES 2 & 3, 13 SEPTEMBER 1973. 


Figure 2, (From ERIM, 1975) 
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Figure 3, Mean radar return as a function of range. The test 

site was divided into 5 strips parallel to the flight line. 
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for each of 1‘he fwo passes (Figures 4 and 5), Passes 2 and 3 were processed separately 
to minimize any biases incorporated intc the data due to variations in aircraft altitude. 
The curves in Figures 4 and 5 are similar in shape (there was no soybeans data for 
Pass 2 beyond range 180) with peaks and valleys occurring at approximately the same 
ranges . 

The radar data were then corrected and means and standard deviations from 
each field were generated. 


4.0 DATA ANALYSIS 

Figures 6 and 7 are scuttergrams of the data for Passes 2 and 3, respectively 
with the two axes being the mean radar return for HH polarization and HV polarization. 
Hyperplanes are drawn on the scattergrams to separate categories . They correspond to 
a discriminant analysis which will be discussed later. The woods give consistently 
high returns, and there is some separability between com and soybeans. The average 
return from all corn and soybeans fields from Passes 2 and 3 were 34.48, 19.46 and 
26.22, 20.96 respectively. Table 1 gives the total number of fields per category for 


each of the two 
PASS 

passes . 
2 

PASS 

3 

Crop Type 

Number of 
Fields 

Crop Type 

Number of 
Fields 

Fallow 

6 

Fallow 

0 

Grains 

5 

Grains 

6 

Pasture 

10 

Pasture 

17 

Woods 

10 

Woods 

16 

Corn 

40 

Corn 

35 

Soybeans 

42 

Soybeans 

42 

Total 

113 

Total 

116 


TABLE 1 . Number of fields per category for 
Poss 2 ond Pass 3 . 
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Flight; Pass 3 
Crop: Soybeans 
Frequency*. L Band 
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Figure 5. Mean radar return for soybeans as a function of range. The test site was divided into 250 strips parallel to the 
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To quonfify ^he measure of separability a linear discriminanf analysis was 
performed on fhe dafa. Fallow and grains v'ere grouped with pasture for two reasons; 

(a) the returns of all three categories were comparable, (b) the grouping increased the 
number of measurements for posture . 

The discriminant analysis uses a regression algorithm which generates a set of 
hyperplanes for separating the training patterns, A total of NC(NC-l)/2 hyperplanes 
a’’e determined, where NC = number of categories . The test patterns are classified on 
th= basis of a majority vote on these hyperplanes [5], For oil the analyses 50% of 
samples were randomly selected to form the training patterns, the remaining samples 
were used as test patterns. Each pattern vector was made up of four measurements— 
like polarization mean and standard deviation and cross polarization mean and 
standard deviation— which shall be represented byHH(M ), HH(a), HV(fi) and HV (cr) 
to facilitate brevity. 

Table 2 shows discriminant analysis results grouping data from both passes for 
four crop types (pasture, woods, corn and soybeans). The probability of correct 
classification is calculated as the number of training + test patterns correctly classified^ 
total number of patterns, HH(M ) is the optimum variable to discriminate crop types if 
only one measurement per pattern vector is used. The optimum combination of two 
variables is HH( p) and HV(fi ) which improves the probability of correct classification 
from 64.6% (for HH ( m )) to 71 .2% (for both HH( M ) and HV( M )) . Table 3 gives the 
contingency tables. The cross polarization component (HV(m)) when added to the like 
polarization return,helps in separating seven fields of woods which were oil classified 
as corn when only the like component was used in the analysis. Some Improvement is 
also noted in the fields of pasture. When woods and pasture were dropped from the 
analysis HH( ) and HH( /i )/HV( ^ ) still gave the best results but both yielded a 
76.1% probability of correct classification. This reinforced the statement that the 
cross polarization return aided only in the separation of woods and pasture from corn 
and soybeans . Adding dimensions to the data did not improve the classification. 

As mentioned previously the data for Pass 2 and Pass 3 were corrected separately 
due to the possibility of system parameters (in particular aircraft height) differing for 
the two passes. A similar analysis os the one mentioned above was conducted separately 
for the two passes. The results are given in Tables 4, 5, 6 and 7. For both passes a 
combination of HH( jJ.) and HV(m) gives the best results end is on improvement of 
about 10% when compared to results using HH(/^) or HV(p) alone . HH(^i) does a good 
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One Dimensional Analysi: 
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67 ouf of 107 fields were correcfly classified. 

Table 2; Discriminanf analysis resulis grouping da^a for bofh passes and using 
four crop types . 















Two Dimensional Analysis (HH(m)/HV(ij) ) 



Table 3: Contingency tables for HH(|i) and HH(^)/HV(^) corresponding 
to Table 2. 
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Table 4: Disrriminanf analysis results for Pass 2 iriing four crop fypes. 


















One Dimensional Analysis (HH(u)) 


Contingency Table for Training Patterns 

Contingency Table for Test Potterns 

True 

True 

Category 

Category 


Woods 


Posfure 0 


Com 


Soybeans 0 



Woods Posture Corn Soybeans 


Assigned Category 


Two Dimensional Analysis (HHIpVHVip)) 


Contingency Table 'or Training Patterns 


Category 


Woods 


Posture 


Com 


Soybeans 



Woods Pasture Com Soybeans 


Assigned Category I 



Woods 


Pasture 


Com 


Soybeans 



1 ° 

1 1 

1 

1 

0 



15 

1 

1 

1 

1 

■ 




Corn Soybec 


Assigned Category 


Contingency Tnble for Test Patterns 


True 

Cotegory 


Woods 


Posture 


Com 


Soybeans I 0 



Assigned Category 


Table 5; Contingency tables for HH(m) and HH(^)/HV(M) corresponding to 
Table 4 (Pass 2). 
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Table 6: Discriminaut analysis resuUs for Pass 3, using four crop types. 


















One Dimensional Analysis (HH(u)) 

Contingency Table for Training Patterns 


0 

0 

"T1 

1 

0 ! 

1 

1 

13 

0 

— 

0 

0 

0 

9 

7 

0 

2 

3 

1 

Woods 

j Pasture 

J 

I Corn 

Soybeans 

1 


Contingency Toble for Test Patterns 


True 

Category 

Woods 

Pasture 

Com 


2 1 

o 

0 ! 

9 1 0 1 

2 

0 ! 8 i 9 

! 1 

0 

1 1 

3 1 4 ' 9 

i 1 

V/oods 

1 Posture 1 Corn 1 Soyoecrs 


Assigr-ed Category 


Two Dimensional Analysis (HH(n)/HV(ui) 

Contingency Table for Training Potterns 


True 

Category 1 
Woods 

Posture 

Com 

Soybeans 



0 

0 

1 

1 

13 1 

0 

0 

1 

1 

1 

1 

0 j 

7 

j 8 

1 

^ i 

3 

1 22 

Woods 

Pasture 

Corn 
— 

1 Soybeans 


Contingency Table for Test Patterns 


Category 


Woods 


Pasture 


C 


Soybeans 



Assigned Category 


Woods 1 Pasture I Corn 1 Soybeor^s 


Assigned Category 


Table 7: Contingency tables for HH(m) and HHIpl/HVlul corresponding 
to Table 6 (Pass 3). 
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job of seporoMng corn and soybeans but confuses woods cs corn and pasture os 
soybeans (for Pass 2). HV ) on the other hand does a better job of separating 
woods ond pasture from corn and soybeans but in terms of total separation it is 
slightly inferior to HH(m ). 

Hypei planes separating categories are drawn on Figures 6 and 7 corresponding 
to the discriminant analysis results given in Tables 5 and 7. 


5.0 CONCLUSIONS 


Analysis of the Huntington County radar data at I -band led to the following 
conclusions: 

1 . The relative radar returns from corn and soybeans at L"band agree 
with the data reported by Ulaby [4] at 4.7 GHz. 

2. It is possible to separate corn, soybeons, woods and pasture with a 
confidence of about 74% if both like and cross polarization returns are 
employed . 

3. If only one polarization is used HH yields good overall rc •♦s (65%) 
and is able to separate cc>n, soybeans and pasture. Howevt woods are 
consistently confused with corn. 

4. The cross polarization component HV is able to differentiate woods 

from the other crop types and if used in conjunction with the like polarization 
component improves the overall confidence of prediction by about 10%. 

Ground based radar data (1-18 GHz) being currently processed and analyzed 
at the University of Kansas were collected over two growing seasons for a variety of 
crop types [6-10] (bare soils, alfalfa, com, soybeans, milo and wheat). The answer 
to the question of what are optimum radar system parameters for the discrimination of 
crops, will be the subject of a future repoi-t. 
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